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Abstract Polarization properties of the photons emitted in the two-photon decay of 
hydrogen-like ions are studied within the framework of the density matrix and second- 
order perturbation theory. In particular, we derive the polarization correlation function 
that gives the probability of the (two-photon) coincidence measurement performed by 
polarization-sensitive detectors. Detailed calculations of this function are performed 
for the 1sii2 ~^ transition in neutral hydrogen as well as Xe^^^ and U^^^ ions. 

The obtained results allow us to understand the influence of relativistic and non-dipole 
effects on the polarization correlations in the bound-bound two-photon transitions in 
heavy ions. 
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1 Introduction 

Studies on the two-photon decay of atomic hydrogen and hydrogen-like ions have by 
now a long tradition both in theory and experiment [l][2]. Originally focused on the 
total transition rates, the interest in these studies has been recently shifted towards 
the energy and angular distributions as well as the 'polarization •properties of the emit- 
ted photons. A series of 7-7 coincidence experiments were performed, for example, to 
measure the polarization correlations in the 2S]^/2 ~^ 1*1/2 decay of \ow-Z ions. When 
compared with theoretical predictions, the experimental results provided an accurate 
test of Bells inequality and hidden variable theories In contrast to \ow-Z sys- 

tems, much less attention has been paid to the polarization properties of the photons 
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Fig. 1 Geometry for the two-photon decay of excited ionic state. 

emitted in the decay of heavy ions. Only in the hght of recent improvements in x-ray 
detector techniques [5], the observation of 7-7 (polarization) correlations in high-Z 
domain becomes more likely within the next few years. Beside the detailed analysis 
of relativistic and quantum electrodynamics (QED) effects, these future coincidence 
experiments might provide an alternative and very promising route for studying the 
parity-violation phenomena in heavy atomic systems or even for testing the symmetry 
violations of Bose particles [E] . 

In this contribution, we apply the density matrix approach in order to analyze the 
spin properties of the photons emitted in the (two-photon) decay of hydrogen-like ions. 
We pay special attention to the so-called polarization correlation function that reflects 
the probability of the 7-7 coincidence measurement performed by two polarization- 
sensitive detectors. The basic relations for such a "polarization-polarization coinci- 
dence" scenario are given in Sections [5] and O In those sections, we discuss how the 
polarization correlation function can be expressed in terms of second-order transition 
amplitude. Even though the evaluation of these amplitudes is a rather complicated task 
due to the summation over the complete spectrum of the system, it can be significantly 
simplified within the non-relativistic dipole approximation for the electron-photon in- 
teraction. Making use of such approximation, we are able then to derive a simple ana- 
lytical expression for the correlation function. In Section |4] we employ this expression 
together with the results from our exact relativistic calculations in order to investi- 
gate polarization correlations in the 1s\i2 ^ ^^1/2 decay of neutral hydrogen as well 
as Xe^3+ and U^^+ ions. Results of our calculations demonstrate that relativistic and 
non-dipole contributions may significantly modify the polarization correlation function 
and that the effect becomes most pronounced for unequal energy sharing between the 
two photons. A brief summary of our results is given finally in Section [S] 



2 Geometry and labeling 

In order to investigate the polarization correlations in the two-photon decay of hydrogen- 
like ions, we shall first agree about the geometry under which the emission of both pho- 
tons is considered. Since, for the decay of unpolarized ionic states, there is no direction 
initially preferred for the overall system, we adopted the z axis along the momentum 
of the "first" photon. Together with the momentum direction of the "second" photon, 
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this axis then defines also the reaction plane (x-z plane). A single polar angle 9, the so- 
called opening angle, is required, therefore, to characterize the emission of the photons 
with respect to each other (cf. Fig. 1). 

As required by Bose-Einstein statistics, the two-photon state has to be symmetric 
upon exchange of the particles. Therefore, it is a priori not possible to address them 
individually. We can safely assume, however, an experimental setup in which two de- 
tectors observe (in coincidence) the photons having certain energies and propagation 
directions. A "click" at these detectors would correspond to the photon's collapse into 
energy and momentum eigenstates [71|S]. A clear identity can be given, therefore, to the 
photons: the first (second) photon is that one detected by the detector A{B) (marked 
gray in Fig. [T]) at a certain energy ijJi{2) and momentum fei(2)- 

For the theoretical analysis below we shall take into account not only the emission 
angles and the energies but also the linear polarization states of the emitted photons. 
In order to observe these states, we assume that both detectors act as linear polarizers 
whose transmission axes are defined in the planes that are perpendicular to the photon 
momenta fci and fc2 and are characterized by the angles xi and X2 with respect to the 
reaction (x-z) plane. In this notation, Xi(2) = 0° denotes the polarization direction of 
the first (second) photon within the reaction plane. 



3 Theory 

Within the framework of the second-order perturbation theory, the analysis of the total 
as well as differential two-photon decay rates is usually traced back to the evaluation 
of the transition amplitude [9lll0]: 



(/|«.u*^e-"-^--|^) (i/|«.u*^e-"-^--'|i) 

E,y — Ei -\- LOl 



+ 



(1) 



where {r\i) = (r) and (r|/) = tpn/jj^ifi'i') are the well-known solutions of the 

Dirac-Coulomb equation for a single electron in the initial and final states, correspond- 
ingly. In this expression, moreover, the transition operator (y.\i\- e^^^'^ describes the 
relativistic electron-photon interaction with u;^ . and being the polarization and the 
wave vector of the j-th photon. 

The evaluation of the amplitude ([T]) is not a simple task owing to the intermediate- 
state summation that includes not only the summation over the discrete part of the 
spectrum but also an integration over the positive as well as the negative-energy con- 
tinuum. A number of methods have been proposed in the past to carry out such a 
summation . In the present work, the second-order transition amplitude ([1} is eval- 
uated by means of the Green's function approach. Since this approach has been widely 
applied over the last years for the analysis of the total two-photon decay rates as well 
as the polarization and angular correlation functions [111112] . it will not be recalled 
here. 



After a brief discussion of the second-order transition amplitude ((T]), we are now 
ready to analyze the polarization properties of the emitted photons. Most naturally. 
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such an analysis can be performed in the framework of the density matrix theory. For 
the decay of an unpolarized initial state \niji) into the level \nfjf), the two-photon 
spin density matrix reads as: 

{ki,Ai,k2,A2|p/lki,A'i,k2,A^) = Yl -^/j'''(^i.^2)A^/l''^*(Ai,A^)(2) 

where Xj = ±1 are the spin projections of the photons onto their propagation directions 
(i.e. the so-called helicity). Instead of this helicity representation, it might be more 
convenient to re-write the density matrix ((2|) in the representation of the vectors Ux 
and My. Such vectors denote the linear polarization of the photons respectively under 
the angles x ~ ^° s^^d x = 90° with respect to the reaction plane (see Fig. [!}. Any 
linear polarization which is nowadays measured in experiments can be expressed in 
terms of these two (basis) vectors 

Ux = (cosxux + sinxuy) = (^e"'^UA=+i + e+'^u^^-i^ , (3) 

by following the standard decomposition of linear polarization vectors in terms of the 
circular polarization ones. 

The density matrix ([2]) still contains complete information on two photons and, 
hence, can be employed to derive their polarization properties. To achieve this, it 
is convenient to define the so-called "detector operator" that projects out all those 
quantum states of the final-state system which lead to a "count" at the detectors. 
Since in our present work we wish to analyze the correlated linear polarization states 
of the photons, the detector operator can be written as: 

P = |kixi k2X2) (kixi k2X2| • (4) 

From this projector operator, by taking the trace over its product with the density 
matrix and applying Eq. ([31, we immediately derive the polarization-polarization 
correlation function: 

^x.x.W-ArTr(Pp,) = ;j^^ ^ 

Ai A2 A2 

which represents the normalized probability of coincidence measurement of two pho- 
tons with well-defined wave vectors ki and k2 and with linearly polarization vectors 
characterized by the angles Xl and X2 with respect to reaction plane. Here the normal- 
ization constant M is chosen in such a way that, for any value of the opening angle 9, 
we get the unity after having summed over the probabilities of the (four) independent 
photons' polarization states \xx), \xy), \yx), \yy). For the sake of brevity, we have in- 
troduced the notation \xx) = |xi = 0°) |x2 = 0°), \xy) = |xi = 0°) |x2 = 90°) and so 
forth. 

Any further evaluation of the polarization correlation ([SJ requires, in general, a 
computation of the fuUy-relativistic transition amplitude ^ . This amplitude accounts 
for the full interaction oc u_\ e*'' '" between the electron and the radiation field and, 
hence, includes the higher non-dipole effects. The non-dipole corrections, however, are 
usually expected to be negligible for low-Z ions. For these ions, it is therefore justified to 
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Zi=0° X,=90° 




9 (deg) 8 (deg) 

Fig. 2 (Color online) Polarization-polarization correlations of piiotons emitted in the 2si/2 ~^ 
two-photon decay. The solid line corresponds to U^^"*", the dot-dashed line corresponds 
to Xe^"^"*" while the dashed line corresponds both to H and to the non-relativistic calculation, 
since the discrepancy between them less is than 10"^"%. Calculations are presented for the 
energy sharing x = 1/16 and for the different angles Xi of the (polarization) transmission axes 
of detectors (cfr. Fig.[T}. 



threat the electron-photon interaction within the non-relativistic dipole approximation 
by setting e''' '^ ~ 1. Within such dipole approximation, a simple analytic expression 
for the function <P can be obtained: 

2E1 1 / \^ 

•^xi .X2 ^ 1 _|_ cos2 g V + COS XI COS X2 cos e j , (6) 

for the particular case of the 2si^2 ^ ^^1/2 two-photon transition. 
4 Results and discussion 

In the present work, we have employed both the non-relativistic electric-dipole approx- 
imation © and the exact relativistic treatment from Eq. (O in order to investigate the 
polarization correlations between the photons emitted in the 2si^2 ~^ ^^1/2 decay. For 
this bound-bound transition, calculations have been performed for neutral hydrogen 
atom H, hydrogen-like xenon Xe^3+ and uranmm 

U91 + 

ions, for various values of the 
energy sharing parameter. This parameter reflects the fraction of energy carried away 
by the "first" photon: x = wi/([.iJi +1^2) and, hence, is defined in the range < x < 1. 
In Fig. 2, for example, the polarization correlation function <?(6) is displayed for the 
parameter x — 1/16. As seen from the figure, for light ions, both the electric dipole and 
the fully relativistic treatments basically coincide and are well described by Eq. ((S)). In 
particular, the polarization correlation function $ almost vanishes if the polarization 
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axes of detectors are perpendicular to each other: xi ~ 0° and X2 = 90° or xi = 90° 
and X2 = 0°. Such a behaviour can be understood if we recall that — within the non- 
relativistic electric dipole approximation — photons are emitted in a pure spin state 
described by the state vector [S]: 

From this expression, it immediately follows that the probability of measuring -in 
coincidence- orthogonal linear polarizations of photons is zero. 

For high-Z hydrogen-like ions, we expect that Eqs. ©-(O might not describe well 
the polarization properties of the emitted photons, owing to relativistic and non-dipole 
effects. As seen from Fig. 2, these effects result in a non-vanishing correlation function 
# for the "perpendicular polarization" measurements, i.e. when X2 ~ Xl ^ 90°. The 
probability of "parallel polarization" measurements consequently decreases of the same 
measure. We moreover notice that, in case of energy sharing x — 1/16 and 6 — 90°, 
events with photon polarizations within the reaction plane are not forbidden, in contrast 
to the non-relativistic dipole approximation, as a consequence of relativistic and non- 
dipole effects. This, in turn, leads to a further reduction of the probability for those 
events with photon polarizations which are perpendicular to the reaction plane. For 
hydrogen-like uranium ions U^^"*", for example, the function ^xi,X2i calculated at 
XI = X2 = 90° and perpendicular photon emission (6 — 90°), decreases from 1 to 
almost 0.85 if the higher multipole terms are taken into account. 

While the relativistic and retardation effects are significant in high-Z domain if 
one of the photons is much more energetic than the second one, they become almost 
negligible if the photons are emitted with nearly the same energy: ui ~ iJ2- As can be 
seen from Fig. 3, for energy sharing x = 0.5, the polarization probabilities obtained 
within non-relativistic electric-dipole and rigorous relativistic approaches differ in fact 
only of about ~ 10^"^, even for the decay of hydrogen- like uranium ions. 

5 Summary and outlook 

In summary, the polarization of the radiation emitted in the two-photon decay of 
hydrogen-like ions has been investigated within the framework of the second-order 
perturbation theory. Special attention has been paid to the correlated spin states of 
two photons as can be measured in coincidence 7—7 experiments. In order to predict the 
outcome of such experiments, the expression for the polarization correlation function 
has been derived within both the exact relativistic theory and the non-relativistic 
electric-dipole approximation. Making use of these two approaches, the photon-photon 
polarization correlations have been calculated for the 2s]^/2 ^ ^^1/2 transition in 
neutral hydrogen as well as hydrogen-like Xe^"^"*" and U^^^ ions. As seen from the 
results obtained, the higher non-dipole terms in the electron-photon interaction may 
affect the correlation function by about 10-20 %; effect that becomes most pronounced 
for the decay of high-.Z ions if the major fraction of the (two-photon) transition energy 
is carried away by a single photon. 

In the present work, we have restricted our theoretical analysis of the two-photon 
decay to the one-electron atomic systems. For high-Z domain, however, the helium-like 
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Zi=0° X,=90° 
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Fig. 3 (Color online) Polarization-polarization correlations of piiotons emitted in the 2si/2 ~^ 
two-photon decay. The solid line corresponds to U^^"*", the dot-dashed line corresponds 
to Xe^"^"*" while the dashed line corresponds both to H and to the non-relativistic calculation, 
since the discrepancy between them less is than 10"^"%. Calculations are presented for the 
energy sharing x = 1/2 and for the different angles Xi of the (polarization) transmission axes 
of detectors (cfr. Fig.[lJ. 



ions are the most suitable candidates for two-photon studies. In the future, therefore, 
we plan to extend our theoretical approach for studying transitions in these two- 
electron species. A first case study on the polarization correlations and spin entangle- 
ment in the lsi/22s]^/2 : J = — >■ ls^^2 '■ J = decay of U^''^ ions is currently under 
way and will be published soon. 
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